Generalized Kadowaki- Woods Relation 
in Heavy Fermion Systems with Orbital Degeneracy 
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We present a theoretical study of the Kadowaki- Woods relation in the orbitally degenerate 
periodic Anderson model. Based on Fermi liquid theory, we derive the generalized Kadowaki- Woods 



relation in the strong coupling limit, A'y 



IQ-y lN{N - 1) [/incm(mol-K/mJ)2], where A is the 



coefRcient of the term in the resistivity, 7 is the T-linear specific heat coefficient, and N is the 
/-orbital degeneracy. This result naturally explains the remarkably smaller value of A'y~'^ in various 
orbitally degenerate (mainly Yb-based) heavy Fermion systems, reported by Tsujii et al., J. Phys. 
Cond. Mat. 15 (2003) 1993. 
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In Fermi liquid (FL) systems, the specific heat C and 
the resistivity p behave as C = 7T and p ^ po + AT^ 
at sufficiently low temperatures. Because 7 oc m* and 
A oc (m*)^ (to* being the effective mass of quasiparti- 
cles) according to the FL theory, the ratio ^7"^ is ex- 
pected to be independent of to*. In fact, many Ce- and 
U-based heavy Fermion (HF) systems follow a universal 
relation, ^7"^ « 1 x 10~^ [/iricm(mol-K/mJ)^], which 
is called the Kadowaki- Woods (KW) relation [1]. More- 
over, this relation also holds in d-electron heavy electron 
systems like LiV2 04 (7 = 350mJ/mole • K^) and A- 15 
compounds (7 ^ 15mJ/mole • K^) [2]. Thus, the KW 
relation has been considered as one of the remarkably ro- 
bust signature of Fermi liquids, irrespective of the value 
of TO*. It can be derived microscopically using the FL 
theory [3], or using the slave-boson method [4]. Note 
that the KW relation holds experimentally even in the 
close vicinity of a magnetic quantum critical point under 
a magnetic field [5], the fact of which is consistent with 
theoretical analyses [6-8]. 

Recently, however, various Fermi liquid systems which 
does not follow the KW relation have been found exper- 
imentally. Especially, Tsujii et al. have revealed that 
^7^^ w 0.4 X 10"^ [/xricm(mol-K/mJ)^] in many Yb- 
based HF systems like YbCu4Ag (7 — 240mJ/mole • K^), 
YbCus^^rAga; (210 - 460mJ/mole-K^), and others 
[2,9,10]. It is about 20~30 times smaller than the con- 
ventional KW ratio, although they are expected to be 
Fermi liquids. Thus, the violation of the KW relation 
should be a very important and fundamental subject on 
the FL theory. The authors of ref. [2] suggest that ma- 
terials with smaller A'y~'^ have almost fully degenerate 
ground states. In fact, the crystalline electric field (CEF) 
in Yb-based HF systems is in general smaller than that 
in Ce-based HF ones because Yb^+-ion is smaller than 
Ce^+-ion, known as lanthanoid contraction. Note that a 
smaller value of ^7^^ is also observed in Pd, Ft, Ni and 
Fe where 7 > 7band [H]- 

In the present work, we revisit the KW relation in 



HF systems based on the FL theory, by taking the /- 
orbital degeneracy into account. By applying the dia- 
grammatic method developed in analyzing the impurity 
Anderson model [12-14], we succeed in deriving a gener- 
alized KW relation, eq. (12), which is valid in the strong 
coupling case where to* 3> TOband- By putting iV = 8 
[N being the orbital degeneracy) which corresponds to 
a free Yb^+-ion, we obtain that ^7"^ ~ 0.36 x lO"*" 
[/ir2cm(mol-K/mJ)^], which is consistent with experi- 
mental observations. This is the first theoretical deriva- 
tion of the KW relation in orbitally degenerate systems. 

As discussed in ref. [2] , the violation of the KW rela- 
tion in Yb-compounds cannot be ascribed to the impurity 
effect nor an accidental singularity of the band structure 
inherent to individual compounds. For this reason, we 
attack this universal issue based on a conventional or- 
bitally degenerate periodic Anderson model (OD-PAM) 
with a single conduction band. Because of the strong L- 
S coupling, the /-electron state for Yb^+-ion with Af^^ 
(Ce+'^-ion with 4/^) is specified by the total angular mo- 
mentum J ~ 7/2 (J — 5/2) and its 2-component Al. 
The degeneracy of /-orbital is A'' = 2 J -I- 1. 

Here, we study the following OD-PAM [15-17]: 
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where cj^^ (/kA/) ^^"^ creation operator of a 
c(conduction)-electron with spin cr = ±i (/-electron 

with M = J,J - l,---,-J), nj^j = fluf^M, Ck 
is the dispersion for c-electrons, and Ei is the local- 
ized /-level energy. Here we assume E{ > p {p 
being the chemical potential) commonly for Ce and 
Yb-based HF systems, by considering the particle-hole 
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transformation for Yb^+-ion. The c- f mixing poten- 
tial matrices are given by VkMo- = K) • <5m,<t for non- 
orbitally-degenerate case (J = 1/2, N = 2), and it is 
given by V^m. = Voy/^^{l /2 + 2Ma)/lYi'^-%ek, ^k) 
for Y h^+ (J = 7/2, iV = 8), and VlcMa = 
2aVb V^V(7/2 - 2Ma)/rYi'^-%6k, ^k) for Ce^+ 
(J = 5/2, iV = 6), where Yj^{9k- (fik) is the spherical har- 
monics. We note that the relation J2m=-j iMiMo-P = Vq 
is satisfied in each case. 

The /-Green function of the present model, GkM m' (e) , 
is given in ref. [16]. It is shown that GkMM'(e) oc 
^i{M-M')<pk^ so it vanishes except for M' = M af- 
ter the summation taken over k. Here we study the 
three-dimensional OD-PAM within the framework of the 
dynamical-mean-field-approximation (DMFA) where the 
d = 00 limit is taken systematically [18,19]. In the 
DMFA, the self-energy I](e) is constructed of local /- 
Green function(s), g{co) = 7^ GkMM(w), which is 
given as [16,17] 

5M = |G'/(a;)+(l- 1)0^0(0;), (2) 

^^('") = E (VG^°M - VoViio + n- eO)"' , (3) 

Gf''{uJ) = {uJ + ^l-Ef-^{uJ))-\ (4) 

where G^^ represents the /-Green function without mix- 
ing with c-clcctrons. A^'o is the number of sites. Note that 
g{u:) is diagonal with respect to M and independent of 
M. Below, we will utilize this fact to analyze the strong 
coupling case. 

First, we discuss the charge susceptibility Xc and the 
specific heat coefficient 7, which are given as 

Xc = -Npf{Q), (5) 

l=~Npf{Q), (6) 

at zero temperature, where \/z = 1 — 9E(e)/9e|c=0j 
1/2^ = 1 - aE(0)/9M, and /(O) = Im.g(-i(5)/7r is the 
/-electron density of states (DOS) per channel at the 
Fermi level; the total /-electron DOS is Np^ (0). 

By shifting the frequencies of every closed loop in E(ti;) 
by w, we obtain the identity [3,14] 

7 = ^ + Eri,Md.M'(0,0)p/(0), (7) 

^ j-M' 

where z,j are site indices, and Ti^M;j,M'{^^^') the full 
four-point vertex within the DMFA; see Fig. 1(a). In 
the strong coupling limit where U is sufficiently large, 7 
will be strongly enhanced whereas Xc is suppressed. This 
means that I/^^m' so 1/^;^ in eq.(7) can be dropped 

in the strong coupling case. 

Here, we introduce a modified renormalization fac- 
tor l/z\oc by dropping Ti^M;i,M' except for j = i in 



eq.(7), which wc call the local approximation. We note 

that by adding the term U Y^^^ ^ImI^Im = U Y.iM ^Im 
to eq.(l) virtually, we can neglect the Pauli principle 
in the Coulomb interaction in constructing diagrams 
[14,16]. This added term can be absorbed by shifting 
Ef in eq.(l) because it is a M-independent constant 
in the paramagnetic uniform state. Then, the identity 
^i,M;i,M' {0;0) = (l-(5MM')rioc(0,0) is easily recognized 
within the DMFA as shown in Fig. 1(b), where rioc(0, 0) 
is the asymmetric local vertex composed of local Green 
functions and U's. We can also check this identity or- 
der by order with respect to U [14]. Taking notice of 
the factor arises from the summation over M' , we find 
that 1/zioc = {N - l)rioc(0,0)p-''(0); see Fig. 1(c). As a 
result, the specific heat coefficient in the local approxi- 
mation, 710C, is given as 

710C = ^N{N - l)rioc(0, 0)/(0)^ (8) 

D 

when 1/zioc 3> 1. We comment that 710c ~ 7 is ex- 
pected in usual paramagnetic heavy Fermion systems 
where magnetic fluctuations are not prominent because 
in such a case the term Ti^M;j,M'{0', 0) with i ^ j, which 
represents the inter-site magnetic correlations, will be 
small [16]. In fact, the universal KW relation A/j'^ ^ 
1 X 10~^/zricm[mol-K/mJ]^ in many Ce- and U-based HF 
systems means the validity of the approximation intro- 
duced in the present work. 

We comment that a similar strong coupling analy- 
sis was performed to derive the Wilson ratio (Wr = 
(x/7)(27T^fc3//i^f})) in the present model, and the rela- 
tion Wr fa 1-1-1/ (Af — 1) is derived on condition that 
the inter-site magnetic-correlations are weak [16]. This 
result will be consistent with a smaller Wr (Wr <, 1) 
in YbCu5_a;Aga;, considering that /Xeff will be slightly 
smaller than 4.54^b (for a free Yb^+ ion) due to a small 
but finite CEF [9]. 

Next, we analyze the imaginary part of the self-energy. 
Its T^-term within the DMFA is give as [3] 

ImSM(O) =ImE(0) = ^^i^ (N - l)rl^{0,0) pf {Of , (9) 

as shown in Fig. 1(d). Using eq. (9), we derive the ex- 
pression for A within the DMFA. According to the Kubo 
formula, the conductivity a is given by 




where /(e) = {e'/^ + 1)"^ and G^(e) = (e -t- /z - Ck - 
Vcf/(e + p - Ef - E(e)))~^ is the Green function for c- 
electrons. Note that vertex corrections for currents arc 
dropped in eq. (10), which is allowed within the DMFA 
[18,19]. 
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The coefficient A = p/T"^ is given by cqs.(lO) and (9). 
Assuming the spherical Fermi surface and using the rela- 
tion (A^/2)/(0) = p^(0)yoV(M-^^/-S(0))2 (p^(0) being 
the DOS for c-electron per spin), we obtain that 



A = 3w'k^^N{N - l)rf„,(0, 0)p/(0)^ 



(11) 



where k-p is the Fermi momentum. The number of elec- 
trons per unit volume in the present model is given by 
n = fcp/37r^. By assuming the free electron model for the 
conduction electrons and reviving h and ks, we obtain 
the following "generalized KW relation" in the strong 
coupling case: 



h 9(37r2)-i/3 



1 



1 X 10-5 



1) 



^N{N-1) 



/iOcm[mol • K/mJ] 



(12) 



where both Fioc and p-^{0) are cancelled out. Here, we 
have used h/e^ = 2.6 x 10^ fl, ks = 1.38 x lO'^^ JK-\ 
and assumed that n~^^^ w a « 1 x 10-^cm (a being the 
lattice spacing). 

According to eq. (12), A'y{^'^ w 1 x lO-^/xOcmfmol • 
K/mJ]^ for N = 2 {J = 1/2), which corresponds to 
the Kramers doublet ground state case due to strong 
CEF. On the other hand, Aj^-J^ w 0.36 x 10-^/incm[mol- 
K/mJ]^ for = 8 (J = 7/2), which corresponds to Yb- 
based HF systems with weak CEF. This result is consis- 
tent with the experiments reported in ref. [2] . 

In the next stage, we study the KW relation in the 
weak coupling region {l/z >, 1) using the second-order- 
perturbation-approximation (SOPA) with respect to U, 
both for J = 1/2 case and J = 7/2 case [17]. In the nu- 
merical calculation, we use the spherical Brillouin zone 
(|k| < tt) for simplicity of the numerical calculation. 
We put ek = -4 + 8(k/7r)^ (the bandwidth being 8), 
Ei = -2.5, Vo = 1.8 and n = 1.15 (n/ = 0.8). Here- 
after, we replace Ej in cq. (1) with Ej + O.OOSek for the 
sake of convenience of numerical calculations. Figure 2 
shows the total /-electron DOS, Npf (0), both for [/ = 
and for i7 = 2 obtained by the SOPA at zero temper- 
ature. The non-interacting DOS for J — 7/2 coincides 
with that for J = 1/2 except for the sharp peak around 
e = Ef ^ 0.6 which is given by G-^"(e) in eq. (2). The 
bottom (top) of the hybridization gap is e = 0.17 (1.84) 
for U = 0. We note again that E[ > /i in the present 
calculation by considering the particle-hole transforma- 
tion for 4/^^-electrons in Yb^+ (J=7/2). In the case of 
U = 2, the DOS around the Fermi level and Ef level are 
renormalized within a smaller energy width due to the 
e-dependence ofReS(e) [17,20]. 

Figure 2 (c) shows the imaginary part of the self-energy 
for U = 2 obtained by the SOPA, which is given by 

lmY,{e-i5)='KU'^{N -I) I du I du' 

Jo J-e+u! 

xpf{oji)pf{J)pf{e-uj + J). (13) 



Thus, Iml]j^7/2(e) coincides with (7/64)ImSj^i/2(e) for 
|e| < \Ef — /u| ~ 0.6 within the SOPA, as shown in Fig. 
2 (c). However, ImSj=7/2(e) takes much larger value for 
|e| >, 0.6 by reflecting the huge weig ht of Gf^{e) aroimd 
Ef. The real part of the self-energy is obtained from 
eq. (13) using the Cauchy integral. Figure 3 shows 
A/7^ as functions of C/^, where A and 7 are obtained 
by the SOPA. [Ah^]j=y^ = (7/64) [A/72] ^^1/2 is re- 
alized within 0(U'^), which holds approximately in the 
weak coupling region where C/^ <, 0.2. We find that 
l/zj=i/2 = 1.25 and l/^j=7/2 = 1.13 for U = 2. Taking 
the result by the SOPA as well as eq. (12) derived by the 
strong coupling analysis, we can naturally estimate that 
the ratio A/j"^ for J = 7/2 is about one order smaller 
than that for J = 1/2 in any intermediate coupling case. 
A perturbation calculation up to i7''-order will be useful 
for a detailed study, which is a future problem. 

Finally, we study the KW relation for the following 
SU(N)-PAM: 



k<T 



kM 



U 



kM 



2 ■-iM''^iM'' (14) 



iM#M' 



where M = J, J - 1, • • • , - J and = 2 J + 1. Ahhough 
this model has been frequently analyzed by slave-boson 
1/A^-expansion method, it is less realistic than eq. (1) in 
that (i) both c- and /-bands have A''-fold degeneracy and 
(ii) the c-f mixing is allowed only for electrons with equal 
M. Apparently, the Green function is diagonal with re- 
spect to M. The local /-Green function g{e) is given by 
G^{e) in eq. (3), instead of eq. (2). 

By performing the same analysis within the DMFA, it 
is shown that eqs.(8) and (9) are also valid in SU(N)- 
PAM in the strong coupling limit. On the other hand, 
the conductivity is given by eq. (10) times A^/2. By us- 
ing the relation pi' (0) = p''{0)V^/{ii - Ef - T,{0) f in 
SU(N)-PAM, we obtain that 



127r^A;p^ 



N- 1 
N 



rL(o,o)p/(o)^ 



(15) 



Considering that fcp = (GTr'^n/N)^/^ in SU(N)-PAM, the 
generalized KW relation in in SU(N)-PAM is given by 



1 X 10 



-5 

J'loc (iV/2)5/3(^ 



1) 



/xOcm[mol • K/mJ]^, (16) 



which is (A^/2)2/'^ times smaller than eq. (12), so it gives 
too small a value for A^ = 8; A'y^'J. ~ 0.14 x lO"'^. We 
note again that SU(N)-PAM is not realistic for larger 
N in that both c- and /-band have A^-fold degeneracy. 
Nonetheless, this result suggests that the degeneracy of 
the conduction band further reduces the value of ^7"^. 
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In summary, we have studied the KW relation in HF 
systems with orbital degeneracy. By analyzing the OD- 
PAM, eq. (1), on the basis of the FL theory, we have 
derived a generalized KW relation in the strong coupling 
limit, eq. (12). The obtained result naturally explains 
the remarkably smaller value of A'y~'^ observed in vari- 
ous Yb-based orbitally degenerate HF systems reported 
in ref. [2]. A numerical analysis using the SOPA was also 
presented. Another generalized KW relation has been de- 
rived based on the SU(N)-PAM, eq. (16), which also tells 
that ^7~^ becomes drastically smaller due to the orbital 
degeneracy. However, the SU(N)-PAM may be less real- 
istic than the OD-PAM for larger N. The present sim- 
plified OD-PAM will be enough to understand a global 
aspect of the KW relation in Ce and Yb-based HF sys- 
tems. It is an important future problem to study the 
effect of small but finite CEF splitting on the value of 
^7"^ by numerical methods. 

Finally, we comment on the pressure dependence of 
A in CeCu2Ge2 [21]: It suggests that the value of A'y~^ 
decreases suddenly when the ground state degeneracy in- 
creases (i.e., AcEF < Tk) for P>15GPa. This interesting 
behavior will be explained within the framework of the 
present study. The generalized K- W relation proposed in 
the present work is confirmed in various HF compounds 
with A'' = 2 ^ 8 [22] . Its importance will increase further 
as various new compounds with orbital degeneracy are 
discovered in future. 

The author is grateful to K. Yamada, T. Saso, D. VoU- 
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and discussions. 
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FIG. 1. (a) a full four-point vertex. (b) a full local 
four-pont vertex, which is given by (1 — &mm')^\oq.- Pioc is 
an asymmetric local vertex, (c) expression for l/^ioc — 1- (d) 
expression for ImE. In (c) and (d), the factor (-N) originates 
from the summation over M' and the Fermion loops. 
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FIG. 2. Total /-electron DOS {Np^ (e)) for (a) (7 = and 
(b) U = 2. A renormalization is recognized in the case of 
U = 2. (c) ImE(e - iS) given by SOPA {U = 2). e = 
corresponds to the Fermi level. 
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FIG. 3. A'y-^ obtained by SOPA as functions of U^. Aj'^ 
for J = 7/2 is about 10 times smaller than that for J = 1/2 
in the weak coupling region. 
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